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Body  of  the  Report 

1 .  Background  and  Content  of  the  Work 

Mixtures  of  polymer  binders  with  new  oxidizers  offer  many  advantages  for  advanced 
propulsions.  Energetic  binders  are  preferable  for  that  propulsions.  The  study  of  combustion 
mechanisms  of  the  energetic  polymer  binders/new  oxidizers  mixtures  is  important,  because 
understanding  of  the  mechanisms  is  useful  for  many  applications  of  these  propulsions. 
Temperature  profiles  of  combustion  waves  were  determined  for  mixtures  with  new  oxidizer  CL- 
20  (HNIW)  and  with  HMX  (70%  of  oxidizer).  Active  binder  PUNE  (plastified  polyurethane 
rubber)  was  used.  The  work  was  performed  at  pressures  0.1-10  MPa  and  sample  temperatures 
from  -100  up  to  +90''C.  Burning  wave  parameters  were  found  by  experiments  and  by  processing 
of  the  experimental  data.  Deep  processing  of  experimental  data  allows  macrokinetics  of 
gasification,  leading  stages,  criteria  of  stable  combustion,  and  sesitivities  of  burning  rate  and 
surface  temperature  to  be  established.  Burn-rate  response  functions  to  oscillatory  pressure  were 
estimated.  Preliminary  results  for  advances  mixtures  (AMMO-BAMO)/HNIW,  20:80,  are 
indicated. 


2.  Approach  and  Object  of  Investigations. 

The  following  propellant  mixtures  were  investigated: 

1.  PUNE/HNIW,  30:70,  p=l,67  g/cm\  sum  formula;  02i(.77H24.u95N2i.62Ci7.697, 

T,d=3171  K(at  10  MPa); 

2.  PUNE/HMX,  30:70,  p=l,7  g/cm^  sum  formula;  O28.51H33.42N21.36C17.566; 

Tod=2717K(at  lOMPa); 

The  used  binder  PUNE  is  a  polyurethane  rubber  (PU),  plastified  by  mixture  of  two  nitroesters 
(NE):  dinitratdietileneglycole  and  dinitrattrietilenglicole.  Relation  of  PU:NE  is  20:80  by  weight. 
Polyurethane  rubber  has  sum  formula:  C48.iH67.67O20.iN2.3.  PUNE  is  an  active  binder,  partly,  due 
to  large  amount  of  nitroesters.  Oxidizer  HNIW  (2,4,6,8,10, 12-hexanitrohexaazaisowurtzitane),  or 
CL-20,  was  a  powder  with  cristal  particles  of  rectangular  shape  of  5x3x3  pm  in  sizes. 
Employment  of  binder  PUNE  with  oxidizer  HNIW  allowed  propellant  samples  having  high 
density,  homogenity  and  high  mechanical  properties  to  be  prepared. 

Mixture  PUNE:HMX,  30:70,  with  high  density,  homogenity  and  mechanical  properties  was 
prepared  to  compare  combustion  mechanisms  of  mixtures  PUNE/HNIW  and  PUNE/HMX. 

3.  Experimental  Techniques. 

3.1.  Conditions  of  experiments  and  microthermocouple  techniques. 

Temperature  profiles  of  the  combustion  waves  and  the  burning  surface  temperatures  weie 
obtained  by  microthermocouple  methods.  Profiles  of  the  combustion  waves  were  obtained  by 
microthermocouples  imbedded  into  solid.  Thermocouples  went  through  combustion  waves  when 
the  waves  propagated  through  the  solid  samples  and  show  temperature  profiles.  The  ribbon  U- 
shaped  thermocouples  made  of  alloys  W+5%ReAV+20%Re  of  2-7  |am  thick  were  imbedded  into 
the  samples.  Cylindrical  samples  of  8  mm  in  diameter  and  3  cm  in  length  were  cat  by  small  thin 
knife  of  U-shape  form,  thermocouple  placed  into  the  created  cross-sections,  which  were  moistend 
by  acetone,  and  the  joined  sections  with  thermocouples  inside  were  protected  by  thin  layer  of 
polymer  glue.  Sample  drying  was  performed  during  6  ours.  Every  sample  had  inside  2-j 
thermocouples  placed  one  above  the  other.  Distances  between  the  junctions  were  2-4  mm.  The 
samples  were  burned  in  a  bomb  of  constant  pressure  in  atmosphere  of  nitrogen  at  pressures  0.1- 
10  MPa  and  at  sample  initial  temperatures  T„=  -100,  +20,  and  +90"C.  In  experiments  at  different 
sample  temperatures,  the  samples  were  placed  into  a  small  thermostat  inside  the  bomb.  Heating  of 
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samples  was  made  by  electrically  heated  wires  and  cooling  was  performed  by  liquid  nitrogen.  The 
sample  temperature  was  controlled  by  thermocouples.  Samples  were  ignited  by  electrically  heated 
wire.  Thermocouple  signals  were  registered  by  amplifier  and  oscillograph.  Burning  rate  was 
measured  by  time  delay  between  the  thermocouple  signals,  by  photoregistrations  of  sample 
combustion  and  by  pressure  increase  during  the  sample  combustion.  Photoregistration  of  sample 
combustion  allows  also  the  type  of  combustion  to  be  established.  Burning  surface  temperatures 
were  measured  by  thermocouples  that  are  being  pressed  to  the  surface  during  sample  combustion 
and  by  establishing  the  locations  of  slope  breaks  on  temperature  profiles  registered  by 
thermocouples  (see  below). 

3.2.  Temperature  profile  measurement  validations.  As  a  rule,  temperature  gradient  close  to  the 
burning  surface  has  a  very  high  value.  It  implies  that  thermocouple  measurements  can  give 
temperature  profiles  with  significant  errors  due  to  thermocouple  heat  inertia.  Because  of  that  it  is 
necessary  to  find  conditions  under  which  thermocouple  measurements  in  combustion  waves  will 
introduce  small  errors.  These  conditions  have  been  found  by  numerical  simulations.  The 
thermocouple  partially  absorbs  the  heat  of  the  thermal  layer  and  decreases  the  temperature  at  the 
point  of  the  measurement.  The  requirement  of  small  temperature  error  (less  than  10%)  of  the 
thermocouple  is  indicated  by  the  following  formula:  h<0.2 

Here:  h  -  thermocouple  thickness  (in  cm),  i-  heat  difflisivity  of  the  solid  (in  cm^/s),  r b  -  linear 
burning  rate  (in  cm/s),  where  /  is  the  thickness  of  heat  layer  of  the  condensed  phase.  The 

thermocouples  have  to  have  U-shape  form.  It  is  necessary  because  of  a  high  difference  between 
heat  conduction  coefficients  of  metallic  thermocouple  wire  and  that  of  solid  or  gas.  Junctions  of 
U-shaped  thermocouples  do  not  experience  large  temperature  decrease.  Modelling  experiments 
and  numerical  simulations  show  that  the  decrease  of  junction  temperature  will  be  small  (<3%)  if 
the  horizontal  part  of  the  U-shaped  thermocouple  is  about  one  hundred  times  more  than 
thermocouple  thickness  h.  Thermal  inertia  of  the  thermocouple  in  the  gas  phase  can  be  taken  into 
consideration  and  eliminated  by  a  correction  procedure.  The  procedure  implies  the  use  of  the 
following  equation: 

dTe,./dx=(rbX„)''-(T  -  Tex); 

Here:  T  -  the  real  temperature  of  gas  in  the  combustion  wave;  Tex  -  the  temperature  registrated  by 
thermocouple;  t„  -  time  response  of  the  thermocouple  in  gas.  Time  response  depends  on  mass 
burning  rate  m  (m=p-rb)  and  temperature  T.  The  temperature  profiles  in  gas  were  corrected  by 
this  equation.  The  theory  of  thermocouple  measurements  in  combustion  waves  of  solids  was 
created  and  confirmed  by  measurements  of  combustion  wave  temperature  profiles  by 
thermocouples  with  sequentially  decreased  thickness  (method  of  “zero  diameter”).  All  the  above 
mentioned  requirements  have  been  met  in  the  investigations.  Different  types  of  metal  wires  for 
thermocouples  were  used  (We, Re  and  Pt,Rh)  to  test  the  catalytic  effect  on  thermocouples.  The 
catalytic  effect  was  not  observed.  The  method  of  burning  surface  temperature  measurement  by 
determining  locations  of  slope  breaks  on  temperature  profiles  registered  by  thermocouples 
(method  of  “slope  break”)  is  based  on  the  existence  of  the  horizontal  delay  of  value  (rbX,,)  on  the 
temperature  profiles  when  thermocouples  go  through  the  burning  surface.  The  delay  is  due  to  the 
change  of  heat  exchange  between  environment  and  thermocouple:  contact  heat  exchange  in  solid 
is  replaced  by  convective  heat  exchange  in  gas. 

4.  Results  of  Measurements  and  Processing 

The  prepared  mixtures  No.  I  (PUNE/HNIW)  and  No. 2  (PUNE/HMX)  were  studied  in  this 
work  detailed  at  sample  temperatures  To=  -  100,  20  and  +90"C.  Uniform  contents  of  the  mixtures 
allow  decreased  scattering  of  all  combustion  characteristics  to  be  obtained.  The  samples  have 
density  p=l,67  g/cm^  for  mixture  No.l  and  p=l,7  g/cm^  for  mixture  No.2.  The  thermodynamic 
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calculated  adiabatic  flame  temperatures  Tad  are  equal  to  3171  K  (at  10  MPa)  and  2717  K  (at  10 
MPa),  correspondingly. 

Mixture  1  and  mixture  2  burn  at  To=  20  and  +90‘’C,  relatively  steady,  without  periodically 
pulsations  of  burning  rate.  Burning  waves  have  stable  propagation  at  these  T„  in  samples  of 
mixture  1  beginning  from  0. 1  MPa,  but  in  samples  of  mixture  2  it  takes  place  only  beginning  from 
0.25  MPa.  At  To=  -  100°C  the  samples  of  mixture  1  burn  steady,  beginning  from  p>  1  MPa  and 
mixture  2  -  beginning  from  p=  1  MPa.  Both  mixtures  burn  relatively  steady  up  to  10  MPa. 

Carbon  residue.  A  significant  carbon  residue  produces  on  the  burning  surface  of  mixture  1  at 
decreased  pressures.  This  residue  approximately  holds  the  shape  of  the  sample  after  going  of  the 
burning  wave  at  pressures  0. 1-1.5  MPa.  At  higher  pressures  the  carbon  residue  has  forms  of 
carbon  flakes  which  accumulale  on  walls  of  the  bomb.  The  mean  density  of  the  residue  pcb  is  low: 
Pcb  =  0.01  -  0.02  g/cm^.  Only  quite  close  to  the  burning  surface  a  thin  layer  of  relatively  dense 
carbon  residue  with  pcb  «  0.1  g/cm^  is  observed  (however,  probably,  this  layer  produces  only 
during  the  sample  extinguishing).  There  is  not  influence  of  To  on  the  carbon  residue  production. 
The  residue,  which  holds  the  shape  of  sample,  is  a  carbon  lattice,  through  which  gaseous 
combustion  products  blow.  These  gases  are  a  good  electrical  insulater  and  because  of  that  the 
thermocouple  signal  don’t  have  any  distortions  when  thermocouple  is  in  the  gas  phase  of  the 
combustion  waves  at  low  pressures. 

Carbon  residue  on  the  burning  surface  of  mixture  2  did  not  observed. 

Temperature  profiles.  Figures  1-12  present  obtained  averaged  temperature  profiles  T(x),  in 
solid  and  gas,  for  mixtures  1  and  2  at  different  pressures  and  sample  temperatures.  Figures  1,  3,  5, 
7,  9,  and  11  present  averaged  distributions  of  temperature  T(x)  along  the  burning  waves  of 
mixture  1  at  0. 1-10  MPa  and  different  To.  Figures  2-12  (even  numbers)  present  the  distributions 
for  mixture  2  at  0.25  -  10  MPa.  The  figures  show  that  both  mixtures  have  temperature  profiles 
T(x)  with  similar  features;  at  0.1-0.25  MPa  the  temperature  profiles  have  only  low-temperature 
(first)  zone  in  the  gas  phase,  at  0.5  MPa  flame  zone  evolve  into  existence  and  at  1.0  MPa  the 
flame  zone  acquires  more  high  temperature.  At  2  MPa  two-zone  sructure  still  exists,  but  at  5-10 
MPa  both  zone  merge  into  one  zone  of  the  gas  phase. 

Temperature  profiles  of  mixture  1  have  significant  temperature  pulsations  in  the  gas  phase 
at  To=  -  lOO^C  and  pressure  1  MPa.  Pressure  increasing  significantly  decreases  amplitude  of  the 
temperature  pulsations  at  To=  -  100“C.  At  To=  20  and  90‘’C  the  temperature  profiles  relatively 
smooth  in  pressure  interval  0.1-10  MPa.  Only  at  pressures  close  to  2  MPa,  increased  temperature 
pulsations  in  the  gas  phase  are  observed  at  the  investigated  To.  A  possible  explanation  of  the 
phenomenon  is  as  follows:  it  takes  place  due  to  carbon  flakes  tearing  off  from  the  burning  surface. 
Indeed,  going  of  that  flakes  close  to  the  thermocouple  junction  have  to  cause  a  local  temperature 
rising  (carbon  lattice  catalize  gas  reactions  and  because  of  that  the  gas  close  to  the  lattice  have  an 
increased  temperature).  At  p>2  MPa  amount  of  flakes  decreases  and,  possibly,  because  of  that 
amplitudes  of  temperature  pulsations  on  temperature  profiles  decrease. 

Temperature  profiles  of  mixture  2  at  all  investigated  pressures  and  initial  temperatures  have  only 
small  irregular  pulsations,  which  are  typical  for  many  propellant  mixtures. 

Processing  of  the  obtained  temperature  profiles  allow  burning  wave  parameters  to  be 
determined.  Tables  1  and  2  present  dependencies  of  the  parameters  on  pressure  at  various  T„  for 
mixture  1.  Tables  3  and  4  present  the  parameters  for  mixture  2. 

Mass  burning  rates  m.  It  can  be  seen  from  Tables  1-4  that  burning  rates  of  mixture  1  are 
significantly  higher  than  the  rates  of  mixture  2.  Indeed  the  burning  rates  of  mixture  1  at  To=  20‘’C 
in  pressure  interval  0.1-10  MPa  comprise  0.22-2.34  g/cm^c,  but  those  of  mixture  2  in  pressure 
interval  0.25-10  MPa  comprise  0.05-  0.87  g/cm^c.  Similar  result  is  observed  at  other  initial 
temperatures.  Indeed,  at  To=  90“C  values  of  m  of  mixture  1  comprise  0.28-2.78  g/cm^c  (0.1-10 
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MPa),  but  for  mixture  2  -  0.08  1,05  g/cm^c  (0.25-10  MPa).  At  T„=  -100“C  m  of  mixture  1 
comprise  0.3-2.12  g/cm^c  (1-10  MPa),  but  for  mixture  2-0.12  0.72  g/cm^c  (1-10  MPa). 

Standard  deviation  of  m  measurements  is  ±5  %. 

Burning  surface  temperatures  Tg.  It  can  be  seen  from  Tables  1-4  that  burning  surface 
temperatures  Tg  of  mixture  1  significantly  lower  than  surface  temperatures  of  mixture  2.  Indeed 
the  surface  temperatures  of  mixture  1  at  To=  20“C  comprise  235-3 18“C  (0.1  -  10  MPa),  but  those 
of  mixture  2  comprise  308-470‘’C  (0.25-10  MPa).  It  can  be  seen  also  that  surface  temperature  of 
mixture  1  increases  slowly  than  that  of  mixture  2  in  the  studied  pressure  interval.  Similar  results 
are  observed  at  other  initial  temperatures.  Indeed,  at  To=  90“C  Tr  of  mixture  1  comprise  243- 
326®C  (0.1-10  MPa),  but  those  of  mixture  2  comprise  320-480“C  (0.25-10  MPa).  At  To=  -100‘’C 
Tg  of  mixture  1  comprise  240-3 12°C  (1-10  MPa),  but  those  of  mixture  2  comprise  330-455“C  (1- 
10  MPa),  Standard  deviation  of  Tg  measurements  is  about  ±5  %. 

Geometric  and  temperature  characteristics  of  temperature  profiles  in  gas. 

As  it  was  mentioned  above,  the  gas  phase  of  the  burning  waves  of  the  studied  mixtures  has  two 
type  of  the  structure:  single-zone  and  two-zone  structures.  The  mixtures  always  have  two-zone 
structure  at  p=0.5-2.0  MPa  and  both  zones  merge  at  p>5.0  MPa.  At  0.1  -0.25  MPa  only  the  first 
flame  exists.  Values  of  mean  temperature  of  the  first  flame  Ti  at  To=  20‘’C  for  mixture  1  are  equal 
to  900-1050‘’C  (0.1-2  MPa)  and  for  mixture  2  -  850-1 100"C  (0.25-2  MPa).  When  initial 
temperature  changes,  values  of  Ti  changes  by  similar  way.  So,  at  To=  OOT  values  of  T|  are  equal 
to  950-1 100°C  (0.1-2  MPa)  and  for  mixture  2  -  1000-1 150“C  (0.25-2  MPa).  At  To=  -100"C 
values  of  Ti  for  mixture  1  are  equal  to  800-900'’C  (1-2  MPa)  and  for  mixture  2  -  950-1000“C  (1-2 
MPa).  Standard  deviation  of  Ti  measurements  is  ±10  %. 

Flame  temperatures  Tf  at  To=  20“C  for  mixture  1  are  equal  to  1600-2860‘’C  and  for  mixture  2  - 
2200-2485‘’C.  When  initial  temperature  changes,  values  of  Tf  also  changes  by  similar  way.  So,  at 
T„=  90‘’C  values  ofTf  for  mixture  I  are  equal  to  1700-2900T  (0.5-10  MPa)  and  for  mixture  2  - 
2250-2500“C  (0.5-10  MPa),  At  To=  -100“C  values  of  Tf  for  mixture  1  are  equal  to  2200-2760“C 
(1-10  MPa)  and  for  mixture  2  -  2150-2330°C  (1-10  MPa), 

Maximal  flame  temperatures  Tf  are  obtained  at  5  MPa  for  both  mixtures  at  all  T„. 

Distances  Li  between  the  burning  surface  and  the  beginning  of  the  flame  zone  decrease  when 
pressure  increases.  In  fact,  Li  is  the  length  of  the  first  flame.  Values  of  Li  at  To=  20‘’C  are  equal  to 
2-  0.4  mm  for  mixture  1  (0.1-2  MPa)  and  to  1.3-0.25  mm  for  mixture  2  (0.25-2  MPa).  Increasing 
To  increases  Li  and  decreasing  To  decreases  Li  for  both  mixtures.  So,  at  To=  90"C  values  of  Li  for 
mixture  1  are  equal  to  2. 5-0.6  mm  (0,1-2  MPa)  aind  for  mixture  2  -  1.5-0. 3  mm  (0.25-2  MPa).  At 
To=  -100"C  values  of  Li  for  mixture  1  are  equal  to  0.4-0.3  mm  (1-2  MPa)  and  for  mixture  2  -  0.4- 
0.2  mm  (1-2  MPa).  The  obtained  data  show  that  the  first  flame  is  an  induction  zone  for  both 
mixtures. 

Distances  L  between  the  burning  surface  and  the  flame  (up  to  0.99  Tf)  also  decrease  when 
pressure  increases.  In  fact,  L  is  the  length  of  the  gas  phase  reaction  zone.  Values  of  L  at  To=  20°C 
are  equal  to  2.5-1. 5  mm  for  mixture  1  (0.5-10  MPa)  and  2. 8-0.7  mm  for  mixture  2  (0.5-10  MPa). 
Increasing  T„  also  increases  L  and  decreasing  T„  decreases  L  for  both  mixtures.  So,  at  T„=  90'’C 
values  ofL  for  mixture  1  are  equal  to  3-1.2  mm  (0,5-10  MPa)  and  for  mixture  2-3-1  mm  (0.5-10 
MPa).  At  To=  -lOOT  values  of  L  for  mixture  1  are  equal  to  3-1.5  mm  (1-10  MPa)  and  for 
mixture  2  -  2. 5-0. 6  mm  (1-10  MPa). 

Standard  deviations  of  Li  and  L  measurements  are  ±15  %. 

Coefficients  of  heat  diffusivity  and  heat  conductivity  in  solid  heat  layer. 

Coefficients  of  solid  heat  diffusivity  were  obtained  by  formula  where  /  is  thickness  of  heat 

layer  of  the  condensed  phase  (see  next  section).  Tables  1-4  show  that  values  of  x  increase  when 
pressure  increases.  Values  of  x  are  independent  on  T„  and  increase  when  mean  solid  heat  layer 
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temperature  increases.  The  values  for  mixture  1  are  equal  to  (l.O-l.l)-lO^  cm^/s  at  0,1  MPa  and 
increase  up  to  (2.7-3. 0)- 10'^  cmVs  at  10  MPa.  For  mixture  2  x  comprise  (0.8-1, 0)- 10'^  cm^/s  at  0.1 
MPa  and  (2.2-2. 5)- 10'^  cmVs  at  10  MPa.  These  values  are  typical  ones  for  x  of  solid  heat  layer  of 
many  propellant  mixtures. 

The  coefficients  of  solid  heat  conduction  were  obtained  by  formula  ;\,s=X-C  p-  It  was  assumed  that 
solid  specific  heat  is  equal  to  C=0.35  cal/g  K  for  mixture  1  and  C=0.33  cal/g  K  for  mixture  2 
(these  values  of  C  were  obtained  by  calculations  of  theory  of  solids).  Tables  1-4  show  that 
increases  with  pressure,  is  independent,  practically,  on  To  and  also  increases  when  mean  solid  heat 
layer  temperature  increases.  Values  of  A-g  comprise  (5.8-17)- lO"^  cal/cm-s-K  for  mixture  1  and 
(4.5-14).  10"^  cal/cm-s-K  for  mixture  2.  These  values  are  also  typical  ones  for  Xs  of  solid  heat  layer 
of  many  propellant  mixtures. 

Characteristic  sizes  in  gas  and  solid. 

Characteristic  thickness  J  and  L  in  solid  (values  of  4  was  detemined  only  for  mixture  2), 
presented  on  Tables  1-4,  were  obtained  by  temperature  distributions  in  the  condensed  phase. 
Experimentally  obtained  thickness  of  heat  layer  in  solid  /  is  a  distance  between  the  burning  surface 
and  a  section  in  solid  with  temperature  T*=(Ts-To)/e  -  To  (e  is  the  base  of  In).  Thickness  of  L  of 
the  partly  melt  layer  (with  melt  HMX)  in  the  condensed  phase  is  a  distance  between  the  burning 
surface  and  a  section  with  HMX’s  melt  temperature  Tn,=280‘’C.  Tables  show  that  thickness  of  / 
and  Im  decrease  when  pressure  increases.  Values  of  /«,  are  less  than  values  of  /. 

Conductive  sizes  &  of  the  gas  phase  were  obtained  by  formula:  S=A,g/mCp;  where  X^  is  coefficient 
of  heat  conduction  in  gas  and  Cp  is  coefficient  of  specific  heat  of  the  gas  phase.  Values  A,g  and  Cp 
were  estimated  by  content  of  combustion  products.  The  obtained  dependencies  of  A,g(Ts)  and 
Cp(Ts)  are  presented  on  Table  5.  Tables  1-4  show  that  values  of  3  decrease  when  pressure 
increases.  Values  of  &  comprise  32-3  pm  for  mixture  1  and  102-6.7  pm  for  mixture  2.  It  can  be 
seen  that  always  S«Li  and  &«L. 

Relative  thickness  of  the  gas-phase  reaction  layer  were  estimated  by  formula  fi=Li/d  at  p<2 
MPa,  and  by  formula  n=L/&  at  p>2  MPa.  It  can  be  seen  from  Tables  1-4  that  inequality  of 
f2»l  always  takes  placee.  It  implies  that  a  wide  reaction  zone  in  the  gas  phase  is  observed  for  all 
combustion  regimes  of  both  mixtures. 

Heat  manifestations:  heat  release  in  solid  and  heat  feedback  from  gas  into  solid. 

The  temperature  profiles  and  burning  surface  temperatures  were  used  for  obtaining  heat 
parameters  of  the  burning  waves.  Heat  flux  from  gas  to  solid  q-m  by  heat  conduction  can  be 
estimated  by  the  following  formula:  q  m=-X.g(T)  (dT/dx)o; 

Here:  (dT/dx)o=q)  is  the  temperature  gradient  in  gas  close  to  the  burning  surface. 

Heat  feedback  q  from  gas  into  solid  by  heat  conduction  can  be  estimated  by  the  following  formula: 

q=-Ag(T)-(p/m: 

Formula  for  heat  release  Q  in  reaction  layer  of  the  condensed  phase  (or  on  the  burning  surface) 
for  mixture  1  has  the  following  form: 

Q~c-(Ts  -To)'q*qr  I 

Formula  for  heat  release  Q  for  mixture  2  is  as  follows: 

Q=c  (Ts  -To)-q-qr+qm: 

where  qm  is  the  heat  of  HMX  melting  in  mixture  2;  temperature  of  HNIW’s  melting  and  heat  of 
HNIW’s  melting  are  not  known  and  because  of  that  formula  for  mixture  1  doesn’t  have  member 
of  qm-  Heat  of  HMX’s  melting  is  equal  to  24  cal/g  (per  gramm  of  HMX)  and  because  of  that 
qm=16.8  cal/g  (per  gramm  of  mixture  2).  Radiant  heat  feedback  qr  from  flame  into  solid  was 
estimated  for  both  mixtures  by  using  content  of  combustion  products  and  real  flame  thick  of  the 
burning  samples  in  our  bomb  of  constant  pressure.  The  heat  radiations  of  CO2,  H2O,  H2,  etc  were 
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taken  into  considerations.  It  was  assumed  that  coefficient  of  heat  radiation  absorbtion  of  the 
burning  surface  was  equal  to  0.5.  Tables  1-4  show  the  obtained  heat  parameters. 

MixturG  1.  It  can  be  seen  from  Tables  1-2  that  mixture  1  has  significant  positive  heat  release  in 
solid.  Values  of  Q  comprise  at  pressures  0.1-10  MPa;  47-90  cal/g  at  To=  20‘’C  and  28-67  cal/g  at 
+90‘’C.  At  pressures  1-10  MPa  and  To=  -lOOT  values  of  Q  are  elevated  and  comprise  91-132 
cal/g.  Parameter  Q  increases  always  with  pressure.  It  can  be  seen  also  that  temperature  gradient  cp 
in  gas  close  to  burning  surface  has  significant  values  and  always  increases  with  pressure. 
Parameter  cp  comprise  in  the  studied  pressure  intervals  values  of  (3-17)-10’  K/cm  at  To=  20“C, 
(3.5-15)  -lO'"  K/cm  at  To=  90“C  and  (5.5-20)  •lO'''  K/cm  at  To=  -100”C.  On  the  contrary,  heat 
feedback  q  from  gas  into  solid  decreases  always  when  pressure  increases.  Values  of  q  are 
comparable  with  values  of  Q  only  at  low  pressures  and  at  To=  +90"C.  Under  other  conditions 
values  of  q  significantly  lowere  than  Q.  Indeed,  Tables  1-2  show  that  q  comprises  in  the  studied 
pressure  intervals  27-15,  25-10  and  31-1 1  cal/g  at  To=  20,  90  and  -lOO'C,  correspondingly. 
Mixture  2.  It  can  be  seen  from  Tables  3-4  that  mixture  2  has  typical  burning  wave  parameters  for 
nitramine  containing  compositions.  Partly,  the  tables  show  that  values  of  Q  are  large  and  values  of 
q  are  relatively  low. 


Table  1:  Burning  wave  parameters  of  mixture  1  (PUNE/HNIW)  at  Tu-ZO^C. 


0.1 

0.5 

1.0 

2.0 

5.0 

8.0 

10.0 

1 

m,  g/cm^s 

0.22 

0.43 

0.62 

0.92 

1.58 

2.06 

2.34 

2 

Ts,  °C 

235 

257 

268 

285 

304 

314 

318 

3 

(p-10-4,  K/cm 

3.0 

5.5 

7.5 

9 

12.5 

15 

17 

4 

q,  cal/g 

27 

26 

24 

18 

16 

15 

16 

5 

qr.  cal/g 

0.8 

0.5 

3.1 

3.3 

4.3 

4.8 

4.8 

6 

Q,  cal/g 

47 

57 

59 

72 

79 

83 

90 

7 

/,  pm 

78 

50 

40 

35 

25 

22 

22 

8 

^(^•10^,  cm^/s 

1.0 

1.3 

1.5 

2.0 

2.4 

2.6 

2.8 

9 

/ts-104, 

cal/cm-s-K 

5.8 

7.5 

8.7 

11.6 

14 

15 

16 

10 

Ti,  °C 

900 

950 

1000 

1050 

*■ 

11 

Tf,  °C 

- 

1600 

2300 

2500 

2860 

2860 

2860 

12 

Li,  mm 

2.0 

i.2 

0.6 

0.4 

- 

- 

“ 

13 

L,  mm 

- 

2.5 

3.0 

2.0 

2.5 

2.0 

1.5 

14 

S,  ^im 

32 

17 

12 

9.5 

5.4 

4.4 

4.0 

15 

Q 

63 

70 

50 

42 

460 

450 

375 

16 

Oo^  kcal/cm^s 

2.4 

8.8 

17.2 

30.6 

70.3 

114 

147 
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Table  2:  Burning  wave  parameters  of  mixture  1  (PUNE/HNIW)  at  To=+90  and  -lOCC 

(+90  OC/-100OC) 


N 

p,  MPa 

0.1 

0.5 

1.0 

2.0 

5.0 

8.0 

10.0 

1 

m,  g/cm^s 

0.28/- 

0.54/- 

0.8/0. 3 

1.2/0.52 

1.98/1.18 

2.5/1.78 

2.78/2.12 

2 

243/- 

265/- 

276/240 

293/260 

312/295 

322/308 

326/313 

3 

(p-lO'"*,  K/cm 

3.5/- 

4.5/- 

6.5/5.5 

8/7 

11/12 

13/17 

15/20 

4 

q,  cal/g 

25/- 

17/- 

16/31 

13/23 

11/18 

10/17 

11/17 

5 

Qr,  cal/g 

0.8/- 

3. 1/3.1 

3.3/3.3 

4. 3/4. 3 

4. 8/4. 8 

4.8/4.8 

6 

Q,  cal/g 

28/- 

44/- 

46/91 

55/100 

62/116 

66/121 

67/132 

7 

/,  nm 

65/- 

43/- 

22/35 

19/25 

18/23 

8 

cm2/s 

1.1/- 

1.4/- 

1.6/1.4 

2.2/1.9 

2.6/2.4 

2,8/2.5 

3.0/2.7 

9 

t>3l- 

8.0/- 

9.3/8.2 

12.7/11 

15/14 

16/14.6 

17/15.8 

cal/cm-sK 

10 

Ti,°C 

950/- 

-/- 

-/- 

-/- 

11 

Tf,°C 

-/- 

1700/- 

2400/2200 

2600/2400 

2900/2760 

2900/2760 

2900/2760 

12 

£ 

E 

2.5/- 

1.8/- 

1.3/0.4 

0.6/0.3 

-/- 

-/- 

-/- 

13 

L,  mm 

-/- 

3/- 

3.5/3 

2.5/2.2 

2.5/2.5 

1. 8/2.2 

1.2/1. 5 

14 

^m 

26/- 

14/- 

9.6/24.8 

6.4/13,8 

4.0/6.8 

3.2/4.8 

3.0/4. 1 

15 

Q 

96/- 

128/- 

135/16 

94/22 

625/368 

560/458 

400/366 

16 

Oo,  kcal/cm^s 

3.6/- 

9.0/- 

19.2/6.2 

35.5/13.5 

80.6/52,4 

120/112 

154/157 

Heat  release  rate  in  gas  close  to  the  burning  surface. 

Heat  release  rates  Oq  in  gas  close  to  the  surface  were  estimated  by  the  use  of  the  obtained  data. 
Oo  shows,  obviously,  the  intensity  of  chemical  reactions  in  gas.  Values  of  Oo  were  obtained  by 
the  following  expression:  Oo=Cp  m-(p; 

Values  of  Cp  were  taken  at  500“C.  Tables  1-4  show  that  values  of  Oq  quickly  increase  with 
pressure.  Mixture  1  has  in  the  studied  pressure  intervals  the  following  values  of  Oq: 

2.4-147  kcal/cm^s  at  To=  20"C,  3.6-154  kcal/cm^s  at  To=  90“C  and  6.2-157  kcal/cm^s  at  To=  - 
100"C.  Mixture  2  has  in  the  studied  pressure  intervals  the  following  values  of  Oo: 

0.28-56  kcal/cm^s  at  To=  20‘’C,  0.64-80  kcal/cm^s  at  To=  90“C  and  1.2-43  kcalWs  at  To=  - 
100°C.  Thus  mixture  1  has  more  intensive  chemical  reactions  in  gas  near  the  burning  surface  than 
mixture  2. 
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Table  3:  Burning  wave  parameters  of  mixture  2  (PUNE/HMX)  at  To-ZO^C. 


No 

p,  MPa 

0.25 

0.5 

1.0 

2.0 

5.0 

8.0 

10.0 

1 

m,  g/cm^s 

0.05 

0.084 

0.16 

0.29 

0.57 

0.76 

0.87 

2 

Ts,  °C 

308 

322 

242 

382 

435 

460 

470 

3 

K/cm 

1.4 

2.1 

3.0 

5.0  . 

10 

14 

16 

4 

q,  cal/g 

50 

45 

34 

30 

32 

32 

34 

5 

Qr.  cal/9 

0.5 

4.4 

3.9 

3.7 

3.3 

2.5 

2.6 

6 

Q,  cal/g 

62 

67 

86 

102 

122 

128 

128 

7 

/,  nm 

270 

220 

140 

90 

60 

50 

50 

8 

Im,  nm 

40 

35 

30 

25 

25 

22 

20 

9 

X’lO^,  cm^/s 

0.8 

1.1 

1.3 

1.6 

2.0 

2.3 

2.5 

10 

ViO^ 

cal/cm-s-K 

4.5 

6.2 

7.3 

9.0 

11 

13 

14 

11 

H 

o 

O 

850 

1000 

1050 

1100 

“ 

12 

H 

o 

O 

- 

2200 

2250 

2400 

2435 

2435 

2435 

13 

Li,  mm 

1.3 

0.9 

0.5 

0.25 

- 

- 

- 

14 

L,  mm 

- 

2.8 

1.8 

2.0 

1.3 

1.0 

0.7 

15 

9,  ^m 

102 

68 

37 

21 

12 

9 

8 

16 

n 

13 

13 

14 

12 

108 

110 

88 

17 

Oq,  kcal/cm^s 

0.28 

0.7 

1.9 

5.8 

22.8 

42.6 

56 
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Table  4:  Burning  wave  parameters  of  mixture  2  (PUNE/HMX)  at  Tq— +90  and  -100"C 

(+90  “CZ-IGOOC) 


N 

p,  MPa 

0.25 

0.5 

1.0 

2.0 

5.0 

8.0 

10.0 

1 

m,  g/cm^s 

0.08/- 

0.14/- 

0.21/0.12 

0.36/0.2 

0.7/0,43 

0.94/0.62 

1.05/0.72 

2 

Ts,°C 

320/- 

340/- 

360/330 

400/360 

450/412 

473/440 

480/455 

3 

(p-lO""^,  K/cm 

2/- 

3/- 

4, 5/2, 5 

7/4 

13/9 

17/13 

19/15 

4 

q.  cal/g 

■  45/- 

39/- 

39/31 

38/32 

33/34 

33/34 

33/36 

5 

Qr,  cal/g 

0,6/- 

3.1/- 

3.5/4 

3.0/4 

2.9/3.5 

2.4/2.7 

2,3/2.7 

6 

Q,  cal/g 

47/- 

58/- 

64/124 

76/133 

99/148 

108/158 

111/161 

7 

/,  urn 

210/- 

140/- 

100/170 

65/120 

45/65 

40/55 

40/50 

8 

Im, 

40/- 

35/- 

30/25 

25/30 

25/22 

25/20 

25/20 

9 

cm^/s 

1.0/- 

1.1/- 

1.2/1.2 

1.4/1.4 

1.8/1.6 

2.2/2.0 

2.5/2.2 

10 

^s-104, 

5.6/- 

6.2/- 

6.7/6.7 

7.8/7.8 

10/9.0 

12/11.2 

14/12.3 

cal/cmsK 

il 

1000/- 

1050/- 

1 100/950 

1150/1000 

7- 

7- 

7- 

12 

Tf,°C 

7- 

2250/- 

2300/2150 

2450/2300 

2500/2330 

2500/2330 

2500/2330 

13 

Li,  mm 

1.5/- 

1.2/- 

0.6/0.4 

0,3/0.2 

7- 

-/- 

-/- 

14 

L,  mm 

7- 

3/- 

2/2.5 

2.5/2.3 

2/1.2 

1. 2/1.0 

1/0.6 

15 

9,  l^m 

71/- 

42/- 

29/45.6 

18/28.6 

10/15 

7.5/10.5 

6.7/9 

16 

o. 

21/- 

29/- 

20/8.8 

17/7.0 

200/80 

160/95 

150/67 

17 

Oo,  kcal/cm^s 

0.64/- 

1.7/- 

3.7/1.2 

10/3.2 

36/15.5 

64/32 

80/43 

Table  5:  Coefficients  of  heat  conduction  Xg  and 
specific  heat  Cp  of  the  gas  phase 


mixtures 

T,  “C 

300 

500 

800 

1000 

1500 

No.  1 

X.g-lO'*, 

cal/cmcK 

1.5 

2.0 

3.0 

3.7 

5.1 

Cn,  cal/g- K 

0.31 

0.37 

0.45 

0.49 

0.62 

No.  2 

X.gl07 

cal/cmc-K 

1.4 

1.6 

2.0 

2.4 

3.8 

Cp,  cal/g- K 

0.38 

0.40 

0.415 

0.42 

0.45 
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5.  Deep  Processing  of  Experimental  Data 

Burning  Rate  Control  Regions. 

The  results  of  measurements  of  q+Pr  and  Q  and  the  investigations  of  the  gas  phase  zone 
behaviour  show  that  burning  rate  of  mixtures  1,  2  are  caused  by  the  heat  release  in  solid  (or  on 
solid  surface)  and,  in  a  smaller  degree,  by  the  heat  feedback  from  the  gas  layer  which  is  close  to 
the  burning  surface.  It  was  mentioned  above  also  that  the  values  of  of  q+qr  decreases  when  the 
burning  rate  increases.  These  facts  show  that  heat  release  in  solid  (or  on  solid  surface)  is  a  main 
factor  of  creating  burning  rate  in  mixtures  1,  2.  Heat  feedback  from  gas  to  solid  is  only  small 
additional  factor  of  the  burning  rate  creating  and  influence  of  this  factor  on  values  of  m  decreases 
when  pressure  increases  (at  pressures  p<10  MPa).  Then,  the  obtained  data,  presented  in  Tables  1- 
5,  allow  an  important  conclusion  to  be  made:  the  high  temperature  region  of  the  burning  waves 
cannot  affect  the  burning  rate.  Indeed,  the  influence  of  the  heat  release  in  the  gas  phase  on  the 
burning  rate  can  be  estimated  by  the  following  formula  obtained  as  a  solution  of  the  heat 
conduction  equation; 

00 

m  q=J  0(x)  exp(-x/fr)-dx; 

0 

This  formula  shows  that  influence  of  heat  release  in  gas  (here  0(x)-is  distribution  of  heat  release 
rate  in  the  gas  phase)  on  m  decreases  very  quickly  when  x>fr.  Tables  1-4  show  that  values  ^  are 
small  and  they  are  significantly  smaller  than  gas  zone  sizes  L  and  Li-  That  implies  that  high 
temperature  region  of  the  burning  waves  cannot  affect  the  burning  rate. 

Thus  the  obtained  results  show  that  the  burning  rate  control  regions  in  the  combustion  waves  for 
all  regimes  of  the  mixtures  combustion,  under  investigated  conditions,  are  the  regions  of  heat 
release  in  solid  just  under  the  burning  surface  (or  immediately  on  the  burning  surface)  and  thin 
low-temperature  gas  layers  close  to  the  burning  surface  (at  x»d).  High  temperature  gas  regions 
cannot  influence  the  burning  rate  in  fact  because  of  a  very  large  heat  resistance  of  the  gas  phase  of 
the  mixtures. 

Macrokinetics  of  Solid  Gasification. 

The  following  equation  connects  burning  rate  of  solid  with  burn-surface  temperature  and 
macrokinetic  characteristics  of  solid  gasification: 

m'=?.sp/Q^-9lTsVE-koQ*- l/Nexp(-E/9tTs) ; 

where  N=l/Tis+(l-%)/Tls-ln(l-%)-(q/Q)  (lii(l-ils)]/%;  ns=Q/Q*;  Q*  is  the  maximum  of 
the  heat  release  in  solid,  E  is  activation  energy  of  limiting  stage  of  the  gasification  process,  Ts  in 
K.  For  studied  combustion  waves  N=l.  The  expression  was  obtained  by  solution  of  system  of 
two  equations  for  the  steady  propagated  burning  wave  (were  solved  the  heat  conduction  equation 
for  solid  phase  and  the  equation  of  diffusion  of  reagents).  The  burning  wave  propagates  due  to 
heat  release  in  solid  Q  and  heat  feedback  from  gas  to  solid  q.  Function  of  the  volumetric  heat 
release  rate  Os  in  solid  was  assumed  as  follows: 

Os(ti,  T)=Q*ko-p-(l-Ti)  exp(-E/9lT): 

Here  r\  is  the  reaction  completeness  (ri=ris  on  the  burning  surface),  ko  is  the  preexponent 
multiplier. 

The  member  exp(-E/9lTs)  plays  the  most  important  role  in  this  formula.  Because  of  that  the 
connection  m  and  Ts  can  be  presented  by  the  following  simplified  expression: 

m= A- exp(- E/29lTs) ; 

This  expression  has  been  designated  “gasification  law”.  The  results  of  measurement,  presented  in 
Tables  1,  2,  allow  the  following  gasification  law  for  mixture  1  to  be  obtained: 
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m=4.0106  exp(-34000/29lTs);  (m  in  g/cm^s,  Ts  in  K) 

Possibly,  value  E=34000  kcal/mole  is  an  activation  energy  of  the  limiting  stage  of  the  process  of 
oxidizer  CL-20  gasification  in  the  reaction  layer  of  mixture  1 . 

The  results  of  measurement,  presented  in  Tables  3,  4,  allow  the  following  gasification  laws  for 
mixture  2  to  be  obtained: 

m=3.8-10^-exp(-25000/29?Ts);  (m  in  g/cm^s,  T  in  K,  at  Ts<470"C), 

The  obtained  value  of  activation  energy  of  the  gasification  process  is  very  close  to  that  of  cyclic 
nitramines  (when  Ts<470“C).  Possibly,  gasification  of  oxidizer  HMX  in  the  reaction  layer  of 
burning  mixture  2  is  the  limiting  stage  of  the  process  of  the  mixture  gasification. 

Pressure  and  Temperature  Sensitivities  of  burning  rate  and  surface  temperature. 

The  following  pressure  and  temperature  sensitivities  for  burning  rate  and  burning  surface 
temperature  were  found  from  the  data  presented  on  Tables  1-4; 

(3=(51nm/5ToVconst,  -  temperature  sensitivity  of  mass  burning  rate; 
r=(5Ts/aToVconsi,  -  temperature  sensitivity  of  burning  surface  temperature, 
v=(51nm/ainp)7b-coi«i,  -  pressure  sensitivity  of  mass  burning  rate; 

|.i=(T»-To)'‘-(aTs/ainp)7b-comt,  -  pressure  sensitivity  of  burning  surface  temperature; 

The  obtained  results  for  mixtures  1  and  2  are  shown  on  Tables  6  for  To=20“C.  It  can  be  seen  that 
burn-rate  pressure  sensitivity  v  for  mixture  1  comprises  0.42-0.58,  and  in  pressure  interval  2-10 
MPa  v=0.58=const.  This  parameter  for  mixture  2  comprises  0.76-0.6.  So,  parameter  v  for 
mixture  1  is  significantly  less  than  that  for  mixture  2.  Burn-rate  temperature  sensitivity  P  of 
mixtures  1  is  equal  to  (0.48-0. 19)- 10'^  1/K  and  that  of  mixture  2  -  (0.3-0. 2)- 10'^  1/K.  So,  burn- 
rate  temperature  sensitivity  of  mixtures  1  is  elevated,  in  comparison  with  that  of  mixture  2. 
Surface  temperature  sensitivity  r  of  mixtures  1  (0.17-0.08)  is  less  than  r  of  mixture  2  (0.17-0.23- 
0.14).  Parameter  \x  of  mixturel  (0.064-0.088-0.064)  also  is  less  than  |a  of  mixture  2  (0.072-0.166- 
0.103). 

Criteria  of  Stable  Combustion.  Criterion  of  stable  combustion  -  Zeldovich’s  criterion 
l<=p.(X^.T^)  .  is  presented  also  on  Table  6.  It  can  be  seen  that  at  pressures  0.25  -  2  MPa  k  can 
exceed  unity.  Novozilov’s  criterion  k*=(k-1)^/(k+l)r  is  necessary  to  calculate  for  regimes  having 
values  of  k>1.  Combustion  is  stable  if  k*<1.  Estimations  of  k*  show  that  mixtures  1  and  2 
always  have  stable  combustion  at  To=20“C.  This  conclusion  confirms  results  obtained  for 
temperature  profile  measurements  -  see  above  n.  4. 

Pressure-Driven  Burn- Rate  Response  Functions. 

The  obtained  results,  which  were  received  under  conditions  of  constant,  nonpulsating  pressures, 
allow  us  to  investigate  burning  rate  behaviour  of  mixtures  1  and  2  under  conditions  of  oscillatoi^ 
pressure.  It  can  be  made  due  to  found  above  sensitivities  for  m  and  Tj.  This  paragraph  is  devoted 
to  these  investigations.  Let  us  assume  that  p=Po‘*‘Pvcoszt7r,  where  pi<<Po  and  uj  is  cyclic 
frequency  of  oscillation  in  1/s.  Then  nondimensional  burning  rate  is  equal  to  V=\+VrCosmt, 
where  V=rJr^o.  V^i=h,7rbo;  r^o  is  a  stationary  burning  rate  and  rb~is  a  pulsating  burning  rate,  which 
takes  place  due  to  pressure  pulsations.  It  is  convenient  to  analyse  burning  rate  response  in 
complex  form.  In  this  case  complex  nondimensional  burning  rate  is  equal  to  V^l+ Kre-'®^  and 
complex  nondimensional  pressure  is  equal  to  7=l+;7re''®^  where  r|i=Pi/po,  and  cy  is 
nondimensional  frequency. 
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The  following  expression  for  burning  rate  response  function  to  oscillatory  pressure  takes  place: 
U=[v+(v-r-p-k)-(z-l)]/[l+r-(z-l)-k-(z-l)/z]; 

Here  z=(l+A/r+4^)/2;  \]=W\It}u  Vi  is  nondimensional  burning  rate  complex  amplitude,  and  o) 
is  nondimensional  frequency  which  is  equal  to  the  cyclic  frequency  m,  in  1/s,  multiplied  by  the 
time  of  solid  heat  layer  relaxation  Values  of  v,  p,  k  and  r  here  are  based  on  the  mean 

pressure  Po  and  because  of  that  data  from  Table  6  can  be  used  for  estimations  of  complex 
functions  of  U.  Thus,  indeed,  we  can  investigate  burning  rate  behaviour  of  mixtures  1,  2  under 
conditions  of  pulsating  pressure. 

Results  of  calculationsof  dimensionless  amplitudes  Re{U} of  burning  rate  pulsations  and  phase  of 
the  pulsations,  Im{U},  as  functions  of  co  are  shown  on  Table  7,  The  used  denominations  are  as 
follows: 

Re{U}n,  is  maximum  of  Re{U}  (in  the  nondimensional  frequency  interval  co  =0-20); 

<y„i  -  nondimensional  frequency,  corresponding  to  this  maximum; 

(U.,2  -  nondimensional  frequency,  corresponding  to  Im{U}  =0; 

Re{U}o  =  Re{U(0)}; 

Re{U}2o=Re{U(20)}; 

Im{Um+}  -  maximum  of  Im{U}  in  interval  cy  =0-<a„2,  i.e.  in  positive  region  of  values  ofImU; 

0)  imi  -  nondimensional  frequency,  corresponding  to  this  maximum; 

Im(Un,-}  -  minimum  of  ImU  (if  this  minimum  takes  place)  in  negative  region  of  ImU; 

CO  in, 2  -  nondimensional  frequency,  corresponding  to  this  minimum; 

G7ni ,  cOra  ,  CTimi  Corresponding  cyclic  frequencies,  in  1/s. 

Estimations  show  that  resonance  frequencies  (y„i,2  are  equal  approximately  to  theoretical  values 
y)„=vT/r.  Typical  response  functions  of  Re{U}  and  Im{U}  for  mixtures  1  and  2  are  shown  on 
Figures  13,  14.  It  can  be  seen  that  response  functions  of  mixtures  1  and  2  are  relatively  close  each 
other.  It  is  natural  that  function  Im{U(<y)},  presenting  the  phase  of  response  functions,  changes 
the  sign  at  the  resonance  frequency.  Functions  of  Re{U((y)}  weakly  depend  on  co,  especially  at 
elevated  p.  The  result  shows  that  resonance  phenomena  in  burning  the  mixtures  are  significantly 
depressed  due  to  viscous  behaviour  of  the  oscillations.  In  fact,  the  pulsating  combustion  of  the 
mixtures  is  an  example  of  a  relatively  smoothed  resonance  phenomenon. 

It  is  necessary  to  indicate  here  same  differences  between  response  functions  of  mixture  1  and  2. 
Table  7  shows  that  maximal  amplitude  of  the  response  Re{U}m  of  mixture  1  in  frequence  interval 
0-20  comprises  1.1-1.5-0.85,  and  that  of  mixture  2  -  2.5-1.4-1.1.  So,  amplitude  reactions  on 
pressure  pulsations  for  mixture  1  are  less  than  that  of  mixture  2.  Table  7  also  shows  that 
differences  (Re{U}„,  -  Re{U}2o),  as  a  rule,  are  less  for  mixture  1,  in  comparison  with  that  of 
mixture  2.  Indeed,  this  value  for  mixture  1  is  0.6-0.7-0.1,  but  for  mixture  2  -  0.8-0.6-0.5.  It  is 
important  to  stress  also  that  differences  (Re{U}2o  -  Re{U}o)  are  positive  for  rnixmre  1  and  they 
are  negative  for  mixture  2.  Both  these  facts  show  that  amplitude  of  response  functions  of  mixture 
1  have  more  viscous  character  than  the  amplitude  of  response  functions  of  mixture  2. 

It  is  a  very  important  result  obtained  here  -  see  Table  7  -  that  cyclic  resonance  frequencies  zsr„i  for 
mixture  1  approximately  an  order  higher  than  cyclic  resonance  frequencies  vJni  fot  mixture  2. 
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Table  6:  Sensitivities  of  burning  rate  and  surface  temperature  of  mixtures  1  and  2  to 
pressure  and  initial  temperature  (for  To=20®C) 


pressure,  MPa 

property 

0.1 

0.25 

0.5 

1.0 

2.0 

5.0 

8.0 

10 

V 

0.42 

0.47 

0.52 

0.58 

0.58 

0.58 

0.58 

r 

0.18 

0.18 

0.17 

0.15 

0.10 

0.08 

0.08 

mixt.l 

0.064 

0.06 

0,08 

0.09 

0.09 

0.08 

0.06 

k 

1.18 

1.19 

1.19 

1.14 

0.85 

0.68 

0.57 

pio-^ 

1/K 

0.55 

0.50 

0.48 

0.43 

0.30 

0.23 

0.19 

V 

0.86 

0.86 

0.86 

0.75 

0.73 

0.62 

0.6 

r 

0.12 

0.14 

0.18 

0.23 

0.21 

0.18 

0.14 

mixt.2 

lA 

0.07 

0.11 

0.17 

0.15 

014 

0.11 

0.10 

k 

1.09 

1.03 

0.98 

1.01 

0.99 

0.97 

0.90 

P•10•^ 

1/K 

0.38 

0.34 

0.30 

0.28 

0.24 

0.22 

0.20 

Table  7:  Response  function  characteristics  for  mixtures  1  and  2  (for  T(,-20"C). 

(presented  as:  mix.  l/mix.2) 


1  pressure,  MPa 

property 

0.1 

0.25 

0.5 

1.0 

2.0 

5.0 

8.0 

10 

Re{U}„, 

1.1/- 

-/2.5 

1. 3/2.0 

1.4/1. 6 

1.5/1. 4 

1. 1/1.4 

0.9/1. 2 

0.8/1. 1 

0)nl 

5.25/- 

-/8.0 

5.5/6.25 

5.5/4.25 

6.0/3.5 

6.5/4.0 

6.8/4.25 

7.3/5.0 

KHz 

0.09/- 

-/0.02 

0.29/0.01 

0.5/0.03 

0.91/0.06 

2.4/0.23 

4.0/0.37 

5.2/0.54 

C0n2 

4.7/- 

-/7.68 

5.2/5.63 

5. 1/3.7 

5. 3/3. 2 

5.4/4.0 

5.65/3.7 

6.25/4.3 

ar„2,  KHz 

0.08/- 

-/0.017 

0.27/0.01 

0.46/0.02 

0.80/0.06 

2.0/0.23 

3.3/0.32 

4.46/0.46 

Re{U}20 

0.5/- 

-/1. 6 

0.6/1.16 

0.66/0.77 

0.79/0.64 

0.83/.66 

0.81/0.58 

0.79/0.64 

Im{U„.} 

0.27/- 

-/0.67 

0.34/0.46 

0.36/0.3 

0.38/0.26 

0.2/0.26 

0.14/0.22 

0,11/0.2 

^  iml 

1.75/- 

-/2.5 

2.0/1.75 

2.0/1.25 

2.0/1.0 

1.5/1.0 

1.25/1.25 

1.25/1.25 

GJiml.KHz 

0.03/- 

-/0.006 

0.1/0.004 

0.18/0.01 

0.3/0.018 

0.56/0.06 

0.74/0.11 

0.9/0.14 

-/- 

-/- 

0.5/0.7 

0.6/0.6 

0.65/0.5 

-/0.5 

-/0.4 

-/0.4 

ini2 

-/- 

-/- 

15/21.5 

15/16 

17/13 

-/14.5 

-/1 6.5 

-/22.5 

^im2y  I'Mz 

-/- 

-/- 

0.78/0.05 

1.3/0. 1 

2.6/0.24 

-/0.8 

-/1. 45 

-/2.43 
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Conclusions 


The  obtained  peculiarities  of  mixture  PUNE/(CL-20),  (30:70),  combustion,  at  To=20"C,  in 
comparison  with  those  of  mixture  PUNE/HMX,  (30:70),  are  as  follows: 

#  a  significant  carbon  residue  on  the  burning  surface  at  pressures  up  to  about  2  MPa  (mixture  2 
has  no  carbon  residue); 

#  significantly  higher  burning  rates,  than  burning  rates  of  mixture  2;  burning  rate  pressure 
sensitivities  are  lower  and  burning  rate  temperature  sensitivities  are  higher  than  those  of  mixtute 
2. 

#  significantly  lower  burning  surface  temperatures,  than  those  for  mixture  2;  surface 
temperature  sensitivities  to  pressure  and  to  initial  temperature  are  also  lower  than  those  for 
mixture  2  (at  p>2  MPa); 

#  heat  release  in  solid  for  mixture  1  is  significantly  lower  than  that  for  mixture  2;  heat  release 
rate  in  gas  close  to  the  burning  surface  for  mixturel  is  significantly  higher  than  that  for  mixture  2; 

#  activation  energy  of  solid  gasifications  in  combustion  waves  is  equal  to  34  Kcal/mole  for 
mixture  1  (it  is  equal  to  24  Kcal/mole  for  mixture  2  in  studied  pressure  interval); 

#  both  systems,  mixtures  1  and  2,  burn  steady,  have  two-zone  structure  of  the  gas  phase  (at 
pressures  up  to  3-5  MPa),  burning  rate  control  for  both  systems  performs  thin  solid  reaction  layer 
on  the  burning  surface  and  thin  low-temperature  gas  layer  close  to  the  surface; 

#  both  systems,  mixtures  1  and  2,  have  resonance  phenomena  in  field  of  pressure  pulsations, 
which  are  significantly  depressed  due  to  viscous  behaviour  of  the  burning  rate  oscillations; 
mixture  1  has  lower  amplitude  of  response  functions,  than  those  of  mixture  2;  absolute  values  of 
cyclic  resonance  frequencies  for  mixture  1  approximately  an  order  higher  than  those  for  mixture 
2. 


6.  Preliminary  results  for  mixture  (BAMO-AMMO)/HNIW,  20:80. 

Mixtures  (BAMO-AMMO)/CL-20,  20:80,  were  prepared  in  the  end  of  this  work. 

Sum  formula:  O23.9H33,05N26.95Ci7.2i;  sample  density  p=l,75  g/cm^ ,  Tad=3200  K  (at  10  MPa). 
Burning  rates  were  obtained  during  preliminary  measurements  at  To=20"C.  Mass  burning  rates 
were  as  follows:  0.6,  2.0  and  5.0  g/cm^s  at  pressures  0.5,  2.0  and  5.0  MPa,  correspondingly. 
Preliminary  measurements  of  temperature  profiles  show,  that  gas  phase  of  the  combustion  waves 
has  two-zone  structure.  The  mesurements  are  continued. 


Recommendations 

In  future  work  must  be  investigate  in  detail  peculiarities  of  combustion  mechanisms  of  mixtures 
on  the  base  of  CL-20,  with  active  binders  GAP,  (AMMO-BAMO),  THF,  and  at  optimal  ratio  of 
oxidizer/binder  (close  to  80:20).  The  investigations  at  different  pressures  and  initial  temperatures 
allow  characteristics  of  the  mixture  combustion  at  stable  pressure  and  at  significant  pressure 
pulsations  to  be  found. 
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temperature  rc  temperature  T'C 


FIGURES 


Figure  1:  Averaged  temperature  profiles  of  mixt.1 
at  0.1  MPa  and  various  To. 


Figure  2:  Averaged  temperature  profiles  of  mixt.2 
at  0.2^  MPa  and  various  To. 
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-0.5  0  0.5  1  1.5  2  2.5  3 


distance  x,  mm 

Figure  3:  Averaged  temperature  profiles  of  mixt.1 
at  0.5  MPa  and  various  To. 


Figure  4:  Averaged  temperature  profiles  of  mixt.2 
at  0.5  MPa  and  various  To. 


17 


temperature  T’C  temperature  T'C 


Figure  5:  Averaged  temperature  profiles  of  mixt.1 
at  1  MPa  and  various  To. 


Figure  6:  Averaged  temperature  profiles  of  mixt.2 
at  1  MPa  and  various  To. 


18 


temperature  T'C  temperature  T'C 


Figure  7:  Averaged  temperature  profiles  of  mixt.1 
at  2  MPa  and  various  To. 


Figure  8:  Averaged  temperature  profiles  of  mixt.2 
at  2  MPa  and  various  To. 
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temperature  TC  temperature  T'C 


distance  x,  mm 

Figure  9:  Averaged  temperature  profiles  of  mixt.1 
at  5  MPa  and  various  To. 


distance  x,  mm 


Figure  1 0:  Averaged  temperature  profiles  of  mixt.2 
at  5  MPa  and  various  To. 
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temperature  TC  temperature  T’C 


Figure  1 1 :  Averaged  temperature  profiles  of  mixt.1 
at  10  MPa  and  various  To. 


Figure  12:  Averaged  temperature  profiles  of  mixt.2 
at  10  MPa  and  various  To. 
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